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We investigate the effects of an external magnetic field in the gluon dynamics of a color supercon- 
ductor with three massless quark flavors. In the framework of gluon mean-field theory at asymptotic 
densities, we show that the long-range component H of the external magnetic field that penetrates 
the CFL phase produces an instability when its strength becomes larger than the Meissner mass 
of the charged gluons. As a consequence, the magnetic field causes the formation of a vortex state 
characterized by the condensation of charged gluons and the creation of magnetic flux tubes. Inside 
the flux tubes the magnetic field is stronger than the applied one. This antiscreening effect is con- 
nected to the anomalous magnetic moment of the gluon field. We suggest how this same mechanism 
could serve to remove the chromomagnetic instabilities existing in gapless color superconductivity. 



PACS numbers: 12.38.Aw, 12.38.-t, 24.85.-|-p 

It is quite plausible [jyj that color superconductivity 
(CS) occurs in the inner regions of neutron stars, hence 
affecting the structure and properties of these compact 
objects. At the same time, it is known that strong mag- 
netic fields exist in the surface of compact stars, reaching 
values of orders 10^'' — 10"'^^ G in the case of magnetars 
[31 • Comparing magnetic and gravitational energies it has 
been shown [3| that the physical upper limit for the total 
neutron star magnetic field is i? ^ 10^* G. If quark stars 
indeed exist and are self-bound rather than gravitational- 
bound this upper limit can be even higher. Therefore, the 
investigation of CS in the presence of a magnetic field is 
of main interest for possible astrophysical applications. 

It is well established by now that at densities high 
enough to neglect the strange quark mass, the ground 
state of three-flavor quark matter corresponds to the 
color-flavor locked (CFL) phase |j| . In this phase, quarks 
form spin-zero Cooper pairs in the color-antitriplet, 
flavor-antitriplet representation. An important feature 
of spin-zero color superconductivity is that although the 
color condensate has non-zero electric charge, the linear 
combination A^ = cos 9 A^ + sin G^of the photon A^ 
and the gluon G® remains massless 0,131 , so it behaves 
as an in-medium electromagnetic field, while the orthog- 



onal combination G^ — 
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cos 9 G,, is massive. 



In the CFL phase the mixing angle 9 is small thus the 
penetrating "electromagnetic" field is mostly formed by 
the original photon with only a small gluon admixture. 

Even though gluons are electrically neutral, in the CFL 
phase they can interact with the modified (also called 
rotated) photon through their Q charges: 
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given in units of e = e cos 9. The Q-charged fields in 
(P) correspond to the combinations G^ = -jk[G'^^ ± «G|^] 
and /^ = -75 [G^ ± iG^]. If the color superconductor is 
penetrated by a sufficiently strong rotated magnetic field, 
the structure and magnitude of the CFL gap are modified 
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iving rise to a new phase called Magnetic CFL (MCFL) 
The MCFL [6j was obtained in the context of an 
effective NJL model in which the interactions between the 
Q-charged gluons and the magnetic field were ignored. In 
this letter, we are interested in extending those previous 
investigations by considering the impact of an external 
magnetic field on the gluon dynamics. We will show that 
at high enough magnetic field there is an inhomogeneous 
condensate of the Q-charged gluons that antiscreens the 
magnetic field due to the anomalous magnetic moment of 
these spin-1 particles. As a consequence, this condensate 
does not give a mass to the Q photon, but instead it 
amplifies the applied Q magnetic field (for the analogous 
behavior of W bosons in a high magnetic field, see Q). 

We will work in the weak coupling region of QCD, 
assuming that the external magnetic field is much smaller 
than the square of the chemical potential, so that the 
effects of the field on the gap structure and magnitude 
can be neglected. In this case, it is safe to assume that 
the superconductor is in the CFL phase. 

We start from the gauge field sector of the mean-field 
effective action in the CFL phase at asymptotic densities. 
This effective action, which includes the one-loop quark 
contribution obtained after integrating out the Nambu- 
Gorkov fields, can be written as 
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where Lg is the gauge fixing Lagrangian density, G" and 
/^i/ are the gluon and electromagnetic field strength ten- 
sors respectively, and 11^^ and II^i, are their correspond- 
ing polarization operators in the CFL phase [SJ . No mix- 
ing between the gauge fields and the Nambu-Goldstone 
mesons is included because this mixing can be eliminated 
using the 't Hooft gauge |2(. 



Since we are interested in the CFL-system response 
to a constant external magnetic field, we can concen- 
trate our attention into the Q-charged gluon part of Q 
in the presence of an external rotated magnetic field H . 
Without loss of generality, the analysis can be done for 
one of the two sets of charged fields, say Gj. For the 
static response, one only needs the leading contribution 
of the polarization operators Q of the Q-charged glu- 
ons in the infrared hmit (po = 0, | p*] — > 0): n°^(a;,y) = 
[rri\)5tj.o5^o + rri\j5fj,i6^i]5{x~y)5°'^, with a,b=4,...,7. The 
Debye (mu) and Meissner {ttim) masses are given by 
ml, = ami, m\j = amy?,, with cr = (21 - 81n2)/18, 
and TOg — g'^[y? /2tt'^). We neglect the possible correc- 
tions produced by the applied field in the infrared masses 
since for fields weak compared to the quark chemical po- 
tential, it is a second order effect. The effective action for 
the fields G^ in the background of the H external field 
becomes 

Kff = Id^^{-liU.r - '^\fL,G- - n,G;|2 



fi ^^ ly 



[(G+)2(G;)2 - iG+G-f] + iG+n^n,G-}, (3) 



where the last term in JSJ comes from the Lg term in |5J), 
taken in the 't Hooft gauge with arbitrary gauge fixing 
parameter A, fl^ = 9^ - ieA^j,, and /^i. = /ff* = H, 
that is, we take the external rotated magnetic field along 
the third spatial direction. The effective action ^ is 
nothing but the characteristic effective action of a spin-1 
charged boson in a magnetic field [iflj. As known, due 
to the anomalous magnetic moment term (ief^i^G'^G'j), 
when the field surpasses a critical value, one of the modes 
of the charged gauge field becomes tachyonic (this is 
the well known "zero-mode problem" found in the pres- 
ence of a magnetic field for Yang-Mills fields I 111 , for 
the W,^ bosons in the electroweak theory |3, U^, and 
even for higher-spin fields in the context of string theory 
[13). The tachyonic mode can be easily found from 
diagonalizing the mass matrix of the field components 
(Gj , G2 ) 
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(4) 



Clearly, above the critical field eHc = rn\j the lowest 
mass mode in (0J becomes tachyonic, with corresponding 
eigenvector of amplitude G in the (1,?) direction for G^ 
(G* in the (1, —i) direction for G"*"). We emphasize that 
the H field producing the instability is the in-medium 
long-range field that does not suffer the Meissner effect 
in the CFL phase. 

Similarly to other spin-1 theories with magnetic in- 
stabilities [7|-[ll|, the solution of the zero-mode problem 
leads to the restructuring of the ground state through the 



formation of a gauge field condensate G, as well as an in- 
duced magnetic field B = V x A due to the backreaction 
of the G condensate on the rotated electromagnetic field. 
The condensate solutions can be found by minimizing 
the Gibbs free energy density Gc — ^ — HB, [T is the 
free energy density), with respect to G and B. Since an 
applied field H > He in the third direction develops an 
instability in the (x, 2/)-plane for the eigenmode G(l, i), it 
is reasonable to make the ansatz Gj~ = — JG^ = G{x,y), 
G^ = Gq = and G^ = (G~)*. Considering this ansatz 
and fixing the gauge parameter to A = 1 in the action 
we obtain for the Gibbs free energy density in the 
G-condensate phase 



Gc = Tno - 2G^n^G - 2{2eB - mlj)\G\^ + 2g^\G\' 



-B^ - HB 



(5) 



where Tno is the system free energy density in the 
normal-CFL phase (G = 0) at zero applied field. 

Using (jSJ the minimum equations for the condensate 
G and induced field B respectively are 

- fl^G - {2eB - mlj)G + 2g^\G\^G = 0, (6) 



2e|Gp -B + H = (7) 

Identifying G with the complex order parameter, Eqs. 
©-0 become analogous to the Ginzburg-Landau equa- 
tions for a conventional superconductor except for the B 
contribution in the second term in (jSJ and the sign of 
the first term in Q. The origin of both terms can be 
traced back to the anomalous magnetic moment term in 
the action of the charged gluons. Notice that because 
of the different sign in the first term of Q, contrary to 
what occurs in conventional superconductivity, the resul- 
tant field B is stronger than the applied field H. Thus, 
when a gluon condensate develops, the magnetic field will 
be antiscreened and the CS will behave as a paramag- 
net. The antiscreening of a magnetic field has been also 
found in the context of the electroweak theory for mag- 
netic fields eH > A/^ 7] . Just as in the electroweak case, 
the antiscreening of the CS is a direct consequence of the 
asymptotic freedom of the underlying theory [7|, [ij] . 

The main goal of this paper will be to investigate the 
qualitative features of the new phase with the gluon con- 
densate at 7? ~ He- A starting point in this direction 
will be to investigate the sign of the condensation energy 



[Gc - Gn]dx 



{-2G^n^G+-{B-Hc) 



-2{2eB - mli)\G\^ + 2g^\G\^}dx (8) 

with G71 = ^no — H^/2 being the Gibbs free energy den- 
sity of the normal-CFL phase in the external rotated 
magnetic field He- Considering that just below He 



the system is in the normal-CFL phase, we should look 
for the possibility of nucleation at He of an inhonioge- 
neous condensate satisfying asymptotic boundary con- 
ditions that correspond to the normal-CFL phase, i.e. 
G{x — > ±00) = 0. Using the Ginzburg-Landau equations 
©-O we can rewrite © as 



a„ 



2 / o 



Hr 



-dx 



(9) 



Because of the hierarchy between the strong {g) and elec- 
tromagnetic (e) interactions, we have that g^ /e^ > 1 and 
hence anc < 0, indicating that the nucleation of an inho- 
mogeneous condensate G is energetically favored. 

To obtain the explicit expression for the condensate 
solution, we can follow Abrikosov's approach 15] to type 
II metal superconductivity for the limit situation when 
the applied field is near the critical value Hc2- Close to 
He the amplitude of the condensate G should be very 
small and we can neglect the nonlinear term in Eq. l|^. 
Similarly, the G^ term in {Tj) can be neglected yielding 
B « He- As a consequence, Eqs. (0-0 decouple and 
G can be found from 
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[d—^Hcxdy~'i7r'-J^x' + -]G = 0, j^x,y (10) 
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where we used the gauge A2 = HqXi and defined $0 



27r/e, and ^^ = l/{2eHc 
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A solution of Eq. (|10|l can be proposed in the form 



G{x,y) = e-'''yyf{x) 
with f{x) satisfying the differential equation 



/"(x)+4^2||(x 



Xkfj^ -j/ 



(11) 



(12) 



and vanishing at a; — > ±00. The parameter x^ is xu = 
ky^o/^nHc- Eqs. l(TU|l and ((T^ are formally identical 
to those found for the Abrikosov's vortex state near the 
critical field Hc2- The parameter ^ in H10|l is playing 
the role of the coherence length, i.e. the characteristic 
length for the variation of G. Increasing the magnetic 
field strength makes ^ to decrease. The solution of (|10|l . 
(O is 



Gk = exp [-iky] exp [- 



2e 



(13) 



where k = ky. From the experience with conventional 
type II superconductivity [la| it is known that the in- 
homogeneous condensate solutions prefer periodic lattice 
domains to minimize the energy. Then, putting on pe- 
riodicity in the y-direction with period Ay = b restricts 
the values of fc to a discrete set k = 27rn/6, n = 1, 2, ... 
This condition implies that we have an infinite set of 



discrete solutions, and because Eqs. H10() and (|12|l are 
linear differential equations, the general solution should 
be given by their superposition G{x, y) — J2 ^nG-n- The 
superposition of all these Gaussian solutions centered at 
different x„ constitutes the vortex state needed to solve 
the instability in the whole space. On the other hand, 
the discrete values of k imply periodicity in x, since the 
Gaussian solutions G„ are located at a;„ = ^"?° =-^^. 

27tHc oHc 

Hence, assuming that all G„ enter with equal weight, 
the periodicity length in the x-direction is Ax — -rj§-- 
Therefore, the magnetic flux through each periodicity 
cell in the vortex lattice is quantized He Ax Ay = <&0i 
with $0 being the flux quantum per unit vortex cell. In 
this semi-qualitative analysis we considered Abrikosov's 
ansatz of a rectangular lattice (i.e. all the coefficients 
C„ being equal), but the lattice configuration should be 
carefully determined from a minimal energy analysis. For 
the rectangular lattice, we see that the area of the unit 
cell is ^ = AxAy = ^q/Hc, so decreasing with H. 

Substituting with Eg. 1)13(1 back into to find a cor- 
rection to the linear solution for the induced field B, one 
easily sees that the induced field, while being homoge- 
neous in the ^-direction, becomes inhomogeneous in the 
(x, 2/)-plane since it depends on the condensate G, which 
has periodicity on that plane. Therefore B forms a flux- 
oid along the z-direction creating a nontrivial topology 
on the {x, y) plane. Notice that the magnetic flux found 
in the present paper is qualitatively different from the 
one found in Ref. [l7|, where the flux was related to the 
massive gauge fleld & . In [l7| the flux was due to the 
existence of nontrivial loops in the order parameter (di- 
quark condensate) space that originated in the transition 
between the unpaired phase and the CFL phase. 

In conclusion, at low H fleld, the CFL phase is an 
insulator, and the H fleld just penetrates through it. 
At sufflciently high H , the condensation of G* is trig- 
gered inducing the formation of a lattice of magnetic 
flux tubes and breaking the translational and remain- 
ing rotational symmetries. Contrary to the situation in 
conventional type-II superconductors, where the applied 
fleld only penetrates through the flux tubes and with a 
smaller strength, the vortex state found in this paper has 
the peculiarity that outside the flux tube the applied fleld 
H totally penetrates the sample, while inside the tubes 
the magnetic fleld becomes larger than H . This antis- 
creening behavior is similar to that of the electroweak 
system at high magnetic fleld 7]. Notice that as the Q 
photons remain massless in the presence of the conden- 
sate G, the C/(l)em symmetry remains unbroken. 

A rough estimate of the critical fleld that produces 
the magnetic instability at the scale of baryon densi- 
ties typical of neutron-star cores (/i ~ 200 — 400Af eT^, 
a,(/i) ~ 1/3) gives He "^ 9.5 x W^G - 3.8 x lO^^G. Al- 
though these are signiflcantly high magnetic flelds, they 
cannot be ruled out as acceptable values for the neutron 



star core. 

We would like to call attention to a possible connection 
between our results and the problem of the instabilities 
[l^|l3 in gapless CS [23 ■ As known, a color supercon- 
ductor can develop chromomagnetic instabilities even in 
the absence of an external magnetic field. These instabil- 
ities may appear after imposing electrical and color neu- 
tralities and (3 equilibrium conditions, and at densities 
where the s quark mass M^ becomes a relevant parame- 
ter. As found first in g2SC [13, and then also in gCFL 
[19|, some charged gluons typically become tachyonic at 
the onset of the gapless phase. From the outcome of this 
paper it results clear that although the imaginary Meiss- 
ner mass of these charged gluons in the gapless phase is 
not triggered by any external magnetic field, the mecha- 
nism we investigated can work here to remove the insta- 
bility of the gapless phase. In this case we could use the 
same effective action (O as a toy model to describe the in- 
stabilities of the charged gluons in gapless CS. The main 
difference will be to consider a negative square Meissncr 
mass from the beginning and to assume that although no 
external magnetic field is present {H = 0), the induced 
magnetic field B can be different from zero (see Eq.(8)) 
since it is nothing but the rotated electromagnetic field 
generated by the gluon condensate. Then the instability 
will be removed through the spontaneous breaking of the 
spatial rotational symmetry SO{3) — > 5*0(2) due to vor- 
tex nucleation and the induction of a rotated magnetic 
field perpendicular to the vortex planes. Notice that the 
new equation for the magnetic field, 2eG^ — B = 0, can 
have a nontrivial solution thanks to the antiscreening ef- 
fect produced by the anomalous magnetic moment contri- 
bution 2eG^. This means that the induction of a S field, 
together with the topological modification of the medium 
that goes with it, can be directly linked to the asymp- 
totic freedom of the theory. It remains to be seen what 
would be the relative field orientation between neigh- 
boring domains that minimizes the system free-energy, 
and whether the generation of G and B alone is enough 
to remove the instability of the neutral gluon modes in 
the gapless phase. If our arguments are corroborated, a 
color superconducting core could provide a new mecha- 
nism to generate and amplify the magnetic fields of com- 
pact stars. The idea of using gluon condensates to solve 
the chromomagnetic instabilities in CS has been also ex- 
plored in Refs. 113. 

It is worth to mention that if the ground state is self- 
consistently found, by minimizing the free-energy with 
respect to all the possible condensates: gluon and diquark 
condensates and the induced rotated magnetic field, the 
inhomogeneous gluon condensate would produce a back- 
reaction into the diquark condensate making it inhomo- 
geneous too and quite likely with the same periodicity of 
the gluon vortex solution. 
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